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O . I. INTRODUCTION. 

<N 

The Casati-Prosen model [1] combines two paradigmatic models of classical and quantum mechanics: a Sinai billiard, 
where the simplest examples of chaotic motion take place and a Young, two slits, experiment, the main example of 
quantum behavior that it "...is impossible absolutely impossible to explain in any classical way" [2]. So we really 
could call this model the " Sinai- Young" experiment. We consider that the complete understanding of this model is 
essential to solve problems like quantum irreversibility, decoherence, and chaos. The model is shown in figure 1 (of 
paper [1]), namely a triangular upper billiard with perfectly reflecting layers, with two slices in its base, on the top of a 
box, the radiating region, with a photographic film in its base and absorbent walls. A quantum state with a gaussian 
packet initial condition bounces in the triangle, and produces two centers of radiation in the two slices from which a 
small amount of probability current leaks from the billiard to the radiating zone. Then when the billiard is perfectly 
triangular and therefore integrable (full-lines in figure 1 of [1]) the interference fringes (full- lines in figure 2 ) appear 
£f) ' in the film and when it is a Sinai billiard and therefore non-integrable (dotted line of fig 1) the first pattern decoheres 
to the (dotted) curve of figure 2. This computer experiment shows how complexity can produce decoherence (without 
an environment or an external noise) and it is explained in paper [1] using a kinematical average. As the subject is 
so important we would like to add another feature to the Casati and Prosen explanation of the phenomenon showing 
■ that the model reaches an equilibrium state where decoherence appears. In doing so we will use our previous results 
Qh! on decoherence [3], mainly paper [4], where local constants of the motion are introduced both at the classical and 
+Jj ■ quantum level allowing to define non-integrable quantum systems and to give a minimal definition of quantum chaos, 
and paper [5], where decoherence times are found. Recently we have shown [6] that our formalism, "Self Induced 
Decoherence" (SID) can be encompassed with the traditional one "Environment Induced Decoherence" (EID) [7], 
combining the advantages of both formalisms. 



> 
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k> ' II. THE PROBLEMS OF PAPER [1]. 

la: 

To make clear our physical point of view let us consider the two main problems to understand chaotic motion in 
terms of quantum mechanics listed in the introduction of paper [1] (see also [8]): 

1. - How is it possible to find chaos in bound systems, with finite number of particles, which have a quasi-periodic 
behavior and therefore a discrete evolution spectrum, if chaotic (e. g. mixing) motion requires a continuous one? 

We consider that the solution can be found in paper [9] where it is shown that, even if a quantum system has a 
discrete evolution spectrum, the motion can be modeled with a continuous spectrum for times much smaller than 
recurrence or Poincare time. For a discrete energy spectrum {a„} this time is 

mm(a„ + i - a v ) 

so if the distances among the eigenvalues are very small tp is extremely large. Then for t <^ tp the typical theorems, 
e.g. the Riemann-Lebesgue theorem, can be used. 

2. - In quantum motions initial errors propagate linearly while in chaotic system this propagation is exponential. 
This contradiction makes quantum chaos impossible. 

We consider that, most likely, this kind of reasonings is done in quantum systems with an integrable classical system 
as classical limit. If this is not the case (as in the systems studied in [4]) it can be demonstrated that the trajectories 
in the classical limit are chaotic and may have positive Lyapunov exponents. So the contradiction is solved. 
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III. THE BASES OF OUR ALTERNATIVE EXPLANATION. 



We will give our alternative explanation based in three results. In this section we only give a sketch of the main 
ideas on these subjects, the complete treatment and figures can be found in the references: 

a.- In paper [4], using the Weyl-Wigner-Moyal isomorphs the definition of classical integrable and non-integrablc 
system is extended to the quantum case. Then the SID formalism is extended to the non-integrable system. For 
N configuration variables these systems, in the classical case have less than N global constants of the motion. But 
according to the Caratheorory-Jacobi lemma [10] they have N constants of the motion locally defined which, via a 
Weyl-Wigner-Moyal isomorphism, allow to define N local Complete System of Commuting Observables that are used 
in the extension of SID. The resulting theory is very similar to the original one. Only an extra index i corresponding 
to the domain D t j >i (that contains the point of the phase space (pi and where local constant of motion are defined) 
must be added in all summations. Then the state of the system p{t) reaches an equilibrium state p*, given by eq. 
(3.23) of [4], defined as a weak limit 



W lim p(t) = p* = / dujp^^^LO, m, m'\ 

,Jo 



(2) 



where to is the eigenvalue of the hamiltonian H (which is considered to be globally defined), {m Xl my) ( j >i in 

our case will be the eigenvalue of the local momentum P^j = (P x , Py)^, and [u),m, m'lfa the cobasis of the eigen 
basis of the CSCO {H^P^.}, namely {\lo, m) t j >i {lo, m!\^ i }. Then the equilibrium final state p» has decohered in the 
energy since only the diagonal terms in lo appear (if not the basis would be (to, lo' , to, to'I^) but not in the remaining 
observables (P x , Py)^, since non-diagonal terms m, ml do appear. Then via a simple diagonalization in the indices 
to, to' we reach to eq. (3.33) 

W lim p{t) = p* = V / dwp(w)^ iP (w,p,p|^ (3) 

t^oo ^ Jo 

where p ( j >i = (p x ,p y )ip i are the eigenvalues of an adequate CSCO {H, O^}. In the correspondent cigenbasis p* is fully 
decohered, since now only diagonal terms (in lo and p) appear. 

b.- The upper triangle will be considered as the Sinai billiard of appendix A of paper [4]. Namely the triangle will 
be complemented by three potential walls in such a way that these potentials U(x,y) (similarly to those of the Sinai 
billiard of appendix A of paper [4]) produce the bounces against the sides of the triangle. We will call D a the interior 
of the triangle (as the D of figure 2 of [4]). It has two independent local constants of the motion: H and P x (or H 
and P y or P^and P y since H — ^y( P x + Py ))- Then we will add three extra domains, each one for each potential 
wall, L>i, D 2 , D4. In the case of the triangle with straight sides the local constants of the motion in the boundaries 
can be deduced by their symmetries. They are: 

• In the horizontal boundary (U(x,y) = U\(y), domain D\) H and P x . 

• In the vertical boundary (U(x,y) = U 2 (x), domain D 2 ) H and P y 

• In the third boundary (U(x,y) — Ui(ax + yb) domain D4) H and a linear combination of P x and P y . 

This is not the case if the triangle has a circular boundary (with radius r — a and angular coordinate 9, where the 
constants of the motion in the third boundary {Ui(x,y) = U(r), domain D4) are H and Pq. 

c- Decoherence times to will be calculated using references [5] and [11]. From reference [5] we know that 

% 

to = - 

7 

where 7 is the distance to the real axis of the pole of the resolvent closer to this axis. These poles for a circular 
symmetric potential can be computed from reference [11] From eqs. (5.1.4) and (5.5.24) of this reference we know 
that the energy is 

E = , and k = — 
2M a 

where M is the mass, being, from eq. (5.5.29), the (3 for the pole closer to the real axis 
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(3o = R - il = U - I ^ ] In - l -ln< ZU ° 



AUq ) \ A 2 ) 2 \ A 2 

where U^{a— ) is the first non vanishing derivative of the potential at the boundary (corresponding to the side of 
the potential) and coefficients Uq and A are given by eqs. (5.5.24) and (5.5.26) of [11]. 
Then 

h 2 R a I a 2Ma 2 
7= ^W< and< ^ Wo (4) 

Below we will use this equation. 



IV. THE TRIANGLE WITH STRAIGHT SIDES. 



Let us first consider the case of the straight triangle and let us take as initial condition in the triangle a pure state 
wave packet \ip) ~ |<^(x,0)} (which of course it is not an eigenstate of the momentum operator P). With this initial 
condition we obtain the solution \tp) ~ \ip(x,t)} in the triangle, that can be written as a matrix 

p(t) = \<p(x,t))(<p(x,t)\ (5) 

We can make some remarks: 

L- If the billiard is considered classical, two initial parallel trajectories remain parallel while they bounce in the 
triangle. Therefore there are neither positive Lyapunov exponent nor chaos. 

ii.- Even if according to eqs. (2) or (3) there will be decoherence in an infinite time, in this case the characteristic 
decoherence time is, in fact, infinite since all the potential walls in this case are straight lines and therefore the radius 
a — > co, then from eq. (4) to — > oo . Therefore p(t) remains bouncing forever in the triangle and does not decohere. 

Let us now consider the lower part under the slit screen.. The direct impact of the packet (5) produces two boundary 
conditions in the two slits:. These two boundary conditions produce two circular-symmetric solutions, <^i(x,i)) and 
|<^2(x,t)), with centers of symmetry in the two slits. Therefore the state in the lower part is \ip(x,t)) = |</?i(x,i)) + 
|i^2(x,t)) and the probability at x is: 

P = {\*){x\)\ v <x,t)){ v <x,t)\ = Pi + Pi + Pint 

where 

pi = Mx,t)| 2 > 0, p 2 = \M*,t)\ 2 > 0, Pint = 2Re(pi(x,t)^(x,t)) (6) 

Of course pi nt ^ and it is the interference term. Let us observe that as \ip(x,t)\ 2 is time invariant p is also time 
invariant as it is verified in [1] figure 4a. 

Moreover if we consider many bounces of the packet instead of just the direct impact, instead of (5) we will have 
a sum with different momenta P. But if the system is integrable this sum will have a finite number of terms (see [1] 
and [12]) and the interference fringes will remain. 



V. THE TRIANGLE WITH A CURVED SIDE. 



Let us now consider the case of the curved triangle. Now 

L- Initial parallel trajectories will lose their parallelism when they collide with the curved side and there will be 
positive Lyapunov exponent and chaos. In fact, Sinai billiards are K-systems. 

ii.- The potential walls are not trivial (i.e. a ^ oo) and therefore the analytic continuation of the resolvent has 
complex poles and the finite decoherence time is given by eq. (4) 1 . 



1 Essentially, as explained, the system has a central domain Do and three potential boundaries domains Di, D2, and D4. But 
each scattering in the potential of D4 can be considered as beginning in the domain Dq = D^ >1 and ending in the out domain 
Dq = D$ 2 . In each scattering the values of the constants of the motion change. As this scattering is repeated again and again 
really Do must be considered as an infinite sequence of D ( p i 
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Then in this case, taking into account the caveat of section II. 1, we can consider that since the system is a K- 
system, it has a continuous spectrum, so using the results reviewed in sections III we can say that the state p(t) inside 
the billiard reaches an equilibrium limit p* given by eq. (2). Now, (ui,p,p'\ ( f )i is a functional (i.e. a distribution or 
kernel) in the continuous variable u> but it is a trivial matrix in the discrete variables p, p'. Nevertheless, based on 
the observation of section II. 1, we can consider u as a discrete variable that has been approximated by a continuous 
one so we can substitute (oj,p 7 p'\ ( p i by \co,p)^ i (u),p'\^, i then the eq. (2) reads 

W lim p{t) = p* = 'Y\p{u)<i >i p\w,p)<i >i (w,p\4, i (7) 

t — >oo * — ' 

iptu 

This is the equilibrium state of the upper part that substitutes the p(t) = |<p(x,i))(<p(x,i)| of eq. (5). Now we must 
obtain the corresponding solution in the lower part, solving the von Neumann equation. But this equation is linear, 
and now the initial conditions are provided not by (5) but by (7), and since in (7) is a linear combination of 
\uj,p) r j >i {uj,p'\ r j >i , to obtain the new p int we must only repeat the same linear combination, e.g. 

Pint = 2 ^] PM« iP WvL P ( x ))#i 

ipu 

where 

<Piwp(x) = <P2wp(x) = (x|w,p} 2 1 

are the solutions in the lower box, centered in the slits 1 and 2 respectively. Now we make the inverse transformation 
the one that allows to go from eq. (2) to eq. (3), i.e. \uj,p) r j >i — U^. \u, m)^ >i where U^. is the unitary transformation 
that diagonalizes (u, m, m'\^ >i so 

ipujmm' 

Y, PMup* [(C^^i wro (x)(^> 2wm '(x))* + (^i wro (x))^> w (x)]^ 

ipujmm' 

Now <^iu;m(x,t) and ip2™'( x '*) are eigenvalues of H and P x , and therefore also of P y , then 2 

But in the slits one of this functions is obtained from the other by a displacement s = (s,0) where s is the distance 
between the slits. Then calling x' = x |s and x" = x + is we have 

^ . m.(x-As) , m'.(x+li) 

Pi„ t = 2 Yl /°("W TO,)**'— «— U = 

ipp' ujmm' 
ipp' ujmm' 

p(ujh iP AU^ i (U^Xe- ii2 ^e i{2 ^ + U^ i (U^Te iSs: =^e- iSsii ^) < i, i 

ipp' ujmm' 

Now we can rephrase what we have said in section II. 1 but now related to the discrete variable m instead of ot v . Now 
the number of terms in the summation is extremely big. In fact, while in the integrable case there will be just a 



2 Really these solutions must be added in order to satisfy the boundary condition of the lower part of the system, but this is 
just another summation that does not modify the final result. 
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few terms (see end of section IV) , now the system is not integrable and the terms may be infinite, since they arrive 
from every direction 3 , so these m are very close. Thus the last J2i P p>mm' can ^ e considered as two integrals in the 
m and in the m' that can be changed in two integrals in the m+ m' and the m— m'. In particular the integrals 

- (m — m' ) - x -(m — m').x 

contain the factors e* } e~ % ^ . So as there is a macroscopic distance from the two slits screen to the 

photographic plate x is macroscopic with respect to h in such a way that we can consider that j- — > oo and we can 
use the Riemann-Lebesgue theorem concluding that 

Pint = j2 p("Wro,(tf;£)*e- ifi ^^ 

ipp ,UJ mm' 

So the interference fringes vanish and there is decoherence in the final equilibrium state, q. e. d. 

VI. CONCLUSION 

1. - We have shown that the Casati-Prosen model reaches an equilibrium state in a finite decoherence time. In this 
final state the interference fringes vanish and we have decoherence. From eq. (4), taking for M the electron mass and 
a = lem we have to ~ Is. 

2. - There is no environment, decoherence is produced by complexity. So the computational result of Casati and 
Prosen cannot be explained by EID. But, we have demonstrated in [6] that a new combined formalism can encompass, 
in a consistent way, EID and SID. In this case SID solves a problem that cannot be solved by EID. The conclusions 
are that EID is a correct theory but it is incomplete and that it can be completed with SID. 

3. - All the reasoning has being done at the quantum level (with some side remarks at the classical level) so we may 
say that the decoherence is produced by quantum chaos. We will try to precise this notion based in SID formalism in 
the near future. 
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